Analysis of the = 4 Hubbard ring using 
counting operators 

Tobias Verhulst* Ben Anthonis, and Jan Naudts^ 

Departement Fysica, Universiteit Antwerpen, 
Groenenborgerlaan 171, 2020 Antwerpen, Belgium 



Abstract 

We prove three theorems about the use of a counting operator to 
study the spectrum of model Hamihonians. We analytically calculate 
the eigenvalues of the Hubbard ring with four lattice positions and apply 
our theorems to describe the observed level crossings. 



1 Introduction 

The Hubbard model is a simplified lattice model used to describe electrons with 
short-range interactions [T] . Its main application is the description of a variety 
of phenomena in solids [HE]- A lot of work has been done on finding exact 
eigenvalues for the Hubbard model and for similar models. In most cases only 
the one dimensional models can be solved analytically. We refer to the one 
dimensional Hubbard model with periodic boundary conditions as the Hubbard 
ring. Most analytic results for the Hubbard ring are found using the Bethe 
ansatz [H O O [3 [3 E] • However, also other methods have been used to analyse 
these models [21 [IHl El [H [I^ . 

One of the interesting aspects of the Hubbard ring is that one observes cross- 
ings between the energy levels which seemingly violate the Wigner-von Neumann 
noncrossing rule [TTl Ul] . One can use the symmetries of the Hamiltonian to 
explain such violations [THIUKin]. Here, we give an explanation which does not 
rely on the integrability of H but rather on that of some reference Hamiltonian 
i?ref- The latter was introduced in [Vl\. This approach is more general in the 
sense that H does not need to be integrable and i?ref can have a much simpler 
spectrum than H. 

We elaborate the method of ^Tf] by introducing a hierarchy of equations. 
Although we will only use this hierarchy to find solutions for the = 4 Hubbard 
ring it should be noted that such a hierarchy always exists. 

Overview of the paper. In Section 2 the Hubbard Hamiltonian is defined and 
some notations are introduced. In Section 3 the hierarchy is introduced and the 
general theorems are proven. Section 4 contains the analytical eigenvalues for 
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the TV = 4 Hubbard ring. We use our theorems to classify these eigenvalues. 
Section 5 gives an overview of the results and conclusions. 

2 Counting operators for the = 4 Hubbard 
ring 

The Hamiltonian of the Hubbard ring is [U [2] 

H = ~Y^tij bl^bj,^ + aYnk,^nk,i (1) 

We consider this model on a lattice with four positions and periodic boundary 
conditions. Furthermore we assume that only nearest neighbour hopping is 
possible and that the band is half filled. This means there are two electrons 
with spin up and two with spin down and the coefficients tij satisfy tij — 
Si.j+i +<5ij_i. Thus the only parameter on which the eigenvalues of ([T|) depend 
is a. 

In [T^, the counting operator, here denoted Mq, is defined as 

Mq = Ynk,ink^l. (2) 
fc 

It counts the number of sites with double occupation. However, the choice of 
counting operator is never unique. It turns out to be advantageous also to 
consider additional counting operators Mi and M2, which differ from Mq by 
adding, respectively, subtracting 1, whenever a doubly occupied site is flanked 
by to singly occupied sites. In formulas. 

Ml = Mo + Va (3) 
M2 = Mo-Va (4) 

where Va is the orthogonal projection onto the space A spanned by the con- 
figurations where two spins on the same site are surrounded by sites with a 
single occupancy. It can easily be verified that [[[i/, Af^], M^], Mj] = [H,Mi] for 
i = 0,1, 2. Then ^17] there exist operators i/iof.i and Ri such that 

H - i7,et,i + R^+ R\ (5) 

with Ri and -fficf,i satisfying 

[F,ef,i,M,] = 

[R„M,] = R, (6) 

For the choices ^ and the corresponding reference Hamiltonians i/rcf,i 
coincide. The extra index i will be omitted in the following, this is, i?rcf = 
Hrcf,i — H-[ci.2- However, i?i does not equal i?2- The method to construct 
explicit expressions for i/rcf en Ri can be found in [17j . It is shown in [17' that, 
in general, the eigenvalues of Mi are of the form = I'n+'m for some constants 
f„. In the iV = 4 Hubbard ring with two spin up and two spin down electrons 
the eigenvalues of the Mi are 0, 1 and 2. 
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3 The hierarchy 



To simplify the notation we will in this section omit the index i. All results 
hold for any of the Mi . In this section we derive a hierarchy of equations which 
is equivalent with the eigenvalue equation. This hierarchy will later be used to 
analyse the spectrum of the A'^ = 4 Hubbard ring. However, the theorems in 
this section are much more general. The only condition for these theorems to 
hold is that the Hamiltonian can be written as ([5]), satisfying ([5]). To simplify 
the calculations we will also assume that the spectrum of AI is of the form 
/i„, = V + m for some v, as is the case for the Hubbard ring. However, the 
theorems can be generalized for models where this is not the case. 
An eigenvector oi H can be decomposed as 

m— 1 

where Mrpm = {y + + m — l)V'm- Here, and V'mmax ^'"^ both non zero, 
Winax can be finite or infinite and mg can be both zero or non zero. Without loss 
of generality we will prove our theorems for mp = 0. Using this decomposition 
of the eigenvector the eigenvalue equation = e'i' can also be decomposed 
into the following hierarchy of equations 

{H,^i-e)ip„, + Rt(;m+i+R''i^m-i=0 (8) 
If TTimax = 1, then ([8|) simplifies to 

iJrofV'l = £1pl (9) 

Now assume that mmax 7^ 1- The first equation of the hierarchy then becomes 

iH,,f - e)^i + R^/J2 = (10) 
If mmax is finite there is also a last equation in the hierarchy which reads 

(i/,ef - £)^^™.„ax + fiV^m„ax-l = (11) 

From ([6]) it is clear that M and RR'^ commute. Thus it is possible to find 
simultaneous eigenstates for AI and RR'^ . The same is true for M and TJrcf- 
However Href does not commute with RR'' . While this means there is no basis 
in which both operators are diagonal it is still possible that there are some 
simultaneous eigenstates. If such states exist they can be used to construct 
special eigenstates of H according to the following Theorem. 

Theorem 1 Let $ be an eigenstate of AI, RR\ (-R^)^ ^'^'^ R with eigenvalues 
/i, and 0. Assume the states $ and 7?^$ are both eigenstates of H^d 
with eigenvalues /i + fc and 11 + a + k + Ak. Then 4' = (1 + /ii?t)$ with 
h = "^^'^ ± ^ ■y/(a + Afc)^ + 4^ are eigenstates of H which have rrimax — 2 
and eigenvalue £ — n + k + hS,. 

Proof 
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Calculating iJ^t gives 

= H{1 + 

= (iJrcf + i?^ + hH.cfR^ + hRR'^) $ 

= (iJrcf + hRR^)^ + (1 + /liJrcf)i?^$ 

= (^ + fc + /iO* + (l + /i(Ai + a + fc + AA:))i?1'<I> (12) 

From this it follows that "if is an eigenstate of H if 

h{fj, + k + h^) = h(p + a + k + Ak) + 1 

This condition is satisfied if h = ± ^ ^^(q! + Akp + 4^. From lEl) follows 

that e = fj, + k + h^. Because -R<E> = = it is straightforward to prove 

that TOmax = 2. 

□ 



In the case of the = 4 Hubbard ring with two spin up electrons and two 
spin down electrons both k and Afc equal zero. However, if an electron is added 
to (or removed from) the model there are states of this form for which Afc does 
not vanish. A similar construction can be done for mmax > 2 states. However, 
the constant h then becomes a solution of higher order equations which means 
there are not always analytical expressions for the eigenvalues. For the N — 4 
Hubbard ring the highest order equation that occurs in these constructions is a 
cubic equation. 

The following theorem about the mmax = 2 case can now be proven (a more 
general form of this theorem, not using the hierarchy, can be found in [T7]). 

Theorem 2 If 'i' is an eigenstate of H for which m^ax — 2, then $ = e^^^^^ 
is also an eigenstate of H for some z £ C, with a different eigenvalue. 

Proof 

Because ^ has Wmax = 2 it can be written as = i/'i + V'2 • Clearly $ can 
then be decomposed as $ = 0i+02 = e~^^tlji + e~^^^tlj2 with R(j)i = R'^(j)2 — 0. 
From Theorem 1 then follows that also {H^ci — aM)(j>i = (i?rcf " aAI)(f>2 = 0. 
Therefore 

H<i> = aM$ + i?Vi +^02 

= vacPi + {va + a)02 + e~^(*^"^^i?t^i + e~^(*^+^'i?V2 

= va(j)i + [va + a)(j)2 — e^{va + a — ^i)4>2 — e^^{va — ^i)4>i 

Hence, $ is an eigenstate of H iff 

(/i — va)e^^ = a + — va — a)e^ (13) 



Thus 



from which z can be calculated. 



va — fj, 
ua — n + a 



(14) 



□ 
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The eigenvalues ei and £2 of two states which are connected in the sense of 
Theorem 2 satisfy 

^ £1+ £2 = a 



£1- £2 = i/a^ + 4^ 



(15) 



with ^ defined as in Theorem 1. 

In the formulation of the third Theorem we use the following notion. The 
hierarchical decomposition ([T]) of the eigenvector 4" is said to be minimal if none 
of the sums 

m' 
m — 1 

with m' < mniax is an eigenvector of H . 

Theorem 3 // $ and ^' are eigenstates of H with the same eigenvalue e. As- 
sume that the hierarchical decomposition of 5* is minimal and that there exists 
a function f such that $ = f{M)^!. Then there exists a constant fi such that 

Proof 

In the case that mmax = 1 nothing has to be proved. Hence assume that 

^max 1 • 

Writing ((TU]) for both 'J and $ = fa{M)'ii gives (to simplify the notation we 
assume 1/ = 0) 

{H,^i - + R^p2 = 0, 

f{^il){H,,t~e)^l + f{^i2)R^2 - 

where /^i, /X2, • • • are the eigenvalues of M. From this pair of equations follows 

if{^ll)-f(P2))m2 = o. 

Note that Rip2 = implies that ipi is an eigenvector of H. But this is not 
possible because of the assumption that the hierarchical decomposition of 5* is 
minimal. Hence, one concludes that /(/ii) = f{l^2)- 

The proof now continues by induction. Assume that /(/ii) = /(/J2) = • • • = 
fiprn') has been proved for some m' < mmax and let us show that this assump- 
tion also holds for m' increased by 1. From ([5]) follows 

(i/ref - £)^rn' + Ri^m' + l + i?Vm'-l = 0, 
f{p7n'){Hrcl - £)lpm + /(/im' + l )i?V'm' + l + f {Uni' -i)R'' Iprn' -1 = 0. 

Using the assumption that f{fJ.m'-i) — fifJ-m') these equations imply that 

(/(A*m' + l) - /(A*m')) Ri^m' + l = 0. 

Now, Rip„i'+i = implies that J2m=i eigenvector of H. This con- 

tradicts the assumption that the hierarchical decomposition of ^' is minimal. 
Hence one concludes that /(/im'+i) = /(Mm')- 

By induction one concludes that all f{pm) are equal. This finishes the proof 
of the Theorem. 



5 



□ 



The Theorem imphes that two eigenvectors $ and of i?, belonging to 
the same multiplet and connected by the operator /(M), cannot have the same 
eigenvalue. This property is expected because of the Wigner-von Neumann non- 
crossing rule. Because M commutes with TJrcf there will be some degeneracies 
in the spectrum of -ffrof which disappear when a perturbation A(i?+i?^) is added 
to i?ref ■ The noncrossing rule states that those eigenvalues which are degener- 
ated for A = do not cross each other for any other value of A. Now, A equals 
one in the Hamiltonian H . Hence, the two eigenvalues of H , which for A = 
are degenerate, must not coincide for A = 1. Note however that a multiplet of 
H must not necessarily originate from a degenerate multiplet of i?ref ■ Hence, 
the result of the Theorem is not an immediate consequence of the Wigner-von 
Neumann noncrossing rule. 

4 The spectrum of H 

Analytical expressions for the eigenvalues as functions of the parameter a have 
been obtained in [T^]. They can also be found using the symbolic software 
package Maple. The eigenvalues are given in Table 1 together with their corre- 
sponding degeneracies. 

Figure 1 shows the eigenvalues of the = 4 Hubbard ring as a function of 
a. It can be seen that there are several crossings between energy levels. The 
crossings at a = are explained by the higher symmetry of the Hamiltonian 
without interaction. However, there are also some crossings for a ^ Q. This is 
similar to the observations of IIJ. 

Using ([7]) the Hilbert space H of eigenfunctions of the Hamiltonian ([T]) can 
be decomposed into subspaces T-Ln,i, where each Mi gives a slightly different 
decomposition. Ri and r\ act as creation and annihilation operators between 
the subspaces mapping TLn^i to T-ln±i.i- In the case of the TV = 4 Hubbard ring, 
the Mi have three eigenvalues (0, 1 and 2) so there are three different subspaces. 
The eigenvalues are then classified by indicating the subspaces Hnfi on which 
which the eigenstates have a non-vanishing orthogonal projection. 

The first three eigenvalues in Table 1 correspond to the states which fall in 
one of the Hn.o- The corresponding eigenstates can easily be constructed. For 
example the eigenvalue 2a corresponds to 

* = <y5l2 - <y5l3 + ^filA. + </'23 - </?24 + ^34 

with the notation ipij = ^b'j ^bl ^bj ^\0) . This eigenvector falls in 7^3.0, the 
eigenvector for which e = falls in Tii.o and the others in 7^2,0- 

The next six eigenvalues, the (1, 2), (2, 2) and (2, 3) pairs, correspond to the 
states that for some n fall within 7i„,o © 'Hn+i.o- The three pairs each form a 
doublet of the form discussed in Theorem 2. Note that there are no eigenstates 
in Hifi Ti-sfi. The reason for this is that if ?/'2 = 0, then R^^i = 0. Thus the 
hierarchy reduces to the single equation ^ and the eigenvalue falls into one of 
the first three classes. 

The remaining eigenvalues correspond to triplet states. The only cases where 
Ml and M2 give different results is for the triplets (1,2,2) and (2,2,3). Using 
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Table 1: The eigenvalues of Hamiltonian ([T]). The functions G±{a) and F±{a) 
are defined by G±{a) = Ta^ ± 36a + 6V-6a4 - 156a2 - 3072 and F±(a) = 



±108a+V-3a6 - lUa* 
with the results of [TT] . 



lOOSa^ — 12288. These eigenvalues are in agreement 



class 



eigenvalue 



degeneracy 



(1) 
(2) 




a 



(3> 


2a 


(1,2) 


\a 




h 16 




\a 


- ya^ - 


h 16 


(2,2) 




a + 2 








a-2 




(2,3) 


\a 


+ §Va2 - 


h 16 



16 



iF+(a)3 + (a2 + i6)F+(a)-^ + 
-\e^F+{a)^ - e-^{a^ + 16)i^+( 



(1,2,3)„ 



(l,2,3)t 



(1,2,2) 



(2,2,3) 



a 



F+{a)^ ~ e^{a^ + 16)^,. (a) 
\F^{a)^ + (a^ + 16)F_(a)-3 +a 

-ie'ri^_(a)* - [a^ + 16)F__(a 
i (G+(a)* + {a^ + 48)G+(a)-* 
~i (e^G+(a)3 +e'T(a2 +48)G+( 
fe'rG+(a)* +eT^(a2 +48)G+( 
i (G_(a)* + (a2 +48)G_(a)-3 +4a) 
-i (e^G_(a)3 +e'3'(a2 +48)G_(a 
-i (e'fG_(a)^ +eT^(a2 +48)G_(a 



a 
a 





the hierarchy derived from Mi one cannot construct the states from the (2,2,3) 
triplet and from M2 one cannot construct the (1,2,2) triplet. 

In Figure 2b it can be seen that within each doublet there are no crossings 
between energy levels. In Figure 2c and 2d it can be seen that this is also 
true for the triplets. This was predicted by Theorem 3. In [T^ it is shown 
that for some models Theorem 2 can be generalized in the sense that from one 
eigenstate $ others can be constructed by defining a function /(M) such that 
^ = f{M)^ is also an eigenstate. Theorem 2 is a special case where f{x) = e^^. 
In [T7] also an example model was given for which there exists triplets of states 
that are connected to each other by such functions and thus the corresponding 
eigenvalues do not cross each other. This is also the case here for the triplets 
in the N = 4 Hubbard model although we do not know the explicit expressions 
for these functions. 
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Figure 1: The complete spectrum of the = 4 Hubbard ring. Both a and e 
are in arbitrary units. Notice the crossings between the levels. Most crossings 
occur for a = but there are also some crossings for other values of a. 

5 Conclusions 

Using the decomposition of the Hamiltonian as described in [17] . we derive a 
hierarchy of equations which is equivalent to the eigenvalue equation. This hi- 
erarchy simplifies the calculation of the eigenstates. Its most important feature, 
however, is that it gives a natural classification of the eigenstates in terms of 
multiplets. 

The first two of our theorems give a method for constructing doublet states 
which are orthogonal to all but two of the subspaces of Ti.. In Theorem 1 a 
method is given to construct the states of such a doublet and Theorem 2 gives 
an expression for the function connecting those states. A construction similar 
to the one of Theorem 1 can also be done for other states, although for the most 
general case the calculations may need to be done numerically. However in the 
A'^ = 4 Hubbard ring there exist only singlets, doublets and triplets which we 
could all identify. Finally, Theorem 3 states that two eigenvalues belonging to 
the same multiplet never cross (for the doublet states this can also easily be 
seen from (fT5|) ). 

In conclusion we use an alternative way to define multiplets without using 
the symmetry of the Hamiltonian. Instead we use operators commuting with a 
reference Hamiltonian i/rcf- Future work also needs to be done on generalizing 
our theorems to triplets and more general multiplets. 
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Figure 2: The eigenvalues of the A'' = 4 Hubbard ring with both a and e in 
arbitrary units, (a) shows the singlets. The solid red line corresponds to (1), the 
dashed blue line to (2) and the black dots to (3). (b) shows the doublet states. 
The pair shown in red dashes is the (1,2) doublet (which falls in Tii ® 'H•l)^ 
the black dotted lines show the (2,2) pair and the solid blue pair is (2,3). In 
(c) (1, 2, 3) a is shown as blue solid lines and (1, 2, 3)^ as red dashes while in (d) 
(1, 2, 2) is shown in blue dashes and (2, 2, 3) in red solid. Notice that within one 
doublet or triplet there are no crossings between different levels. 
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